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ABSTRACT 



■ A spectroscopic technique employing an iodine absorption cell (I2) to super- 

impose a reference spectrum onto a stellar spectrum is currently the most widely 
CSj ' adopted approach to obtain precision radial velocities of solar-type stars. It 

^ 1 has been used to detect ~80 extrasolar planets out of ~130 known. Yet in its 

t> ■ original version, it only allows us to measure precise radial velocities of single 



stars. In this paper, we present a novel method employing an I2 absorption cell 



OO 

that enables us to accurately determine radial velocities of both components of 
double-lined binaries. Our preliminary results based on the data from the Keck I 
telescope and HIRES spectrograph demonstrate that 20-30 ms" 1 radial veloc- 
ity precision can be routinely obtained for "early" type binaries (F3-F8). For 
Qh 1 later type binaries, the precision reaches ~10 ms -1 . We discuss applications of 

the technique to stellar astronomy and searches for extrasolar planets in binary 
systems. In particular, we combine the interferometric data collected with the 
Palomar Testbed Interferometer with our preliminary precision velocities of the 
spectroscopic double-lined binary HD 4676 to demonstrate that with such a com- 
^ ■ bination one can routinely obtain masses of the binary components accurate at 

■ least at the level of 1.0%. 

Subject headings: binaries: spectroscopic — stars: fundamental parameters - 
stars: individual (HD 4676, HD 109358, HD 206901, HD 209458, HD 282975) - 
techniques: radial velocities 



1. Introduction 

The idea of passing starlight through an absorption medium to superimpose a set of 
reference lines and subsequently precisely measure radial velocities was first proposed by 
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Griffin & Griffin (1973). At that time, Griffin & Griffin (1973) used the atmosphere and 
the telluric H 2 and 2 lines. In order to overcome systematic errors in velocities caused 
by, among others, the atmospheric pressure changes and hence shifts in the position of the 
telluric lines, Campbell & Walker (1979) introduced for the first time an absorption cell 
(containing a toxic hydrogen fluoride gas). This technique was subsequently used to carry 
out a pioneering precision radial velocity survey of 16 stars over the period of six years (1981- 
1987; Campbell, Walker, & Yang 1988). It is conceivable that Campbell, Walker, & Yang 
would have detected the first extrasolar planet if they targeted a larger number of stars as 
their average velocity precision was ~13 ms _1 . 

Marcy & Butler (1992) took the absorption cell technique to a new level by introducing 
an iodine (I 2 ) absorption cell with the modeling of a spectrograph point spread function 
(PSF). The gaseous iodine, besides being much less toxic, has a strong line absorption 
coefficient (and hence requires a path length of only a few centimeters) and offers a dense 
forest of absorption lines over the wavelength range of 500-630 nm. In their approach, the 
Doppler shift of a star spectrum AA is determined by solving the following equation (Marcy 
& Butler 1992) 

Iobs(X) = + AA S )T /2 (A + AA /2 )] ® PSF (1) 

where AA S is the shift of the star spectrum, AA/ 2 is the shift of the iodine transmission 
function Tj 2 , <g) represents a convolution and PSF a spectrograph PSF. The parameters 
AA S , AA/ 2 as well as parameters describing the PSF are determined by performing a least- 
squares fit to the observed (through the iodine cell) spectrum I obs . To this end, one also 
needs a high SNR star spectrum taken without the cell I s which serves as a template for all 
the spectra observed through the cell and the I 2 transmission function Tj 2 obtained with the 
Fourier Transform Spectrometer at the Kitt Peak National Observatory. The Doppler shift 
of a star spectrum is then given by AA = AA S — AA/ 2 . 

To date, the iodine technique has been successfully used to detect ~80 out of ~130 
known extrasolar planets (Schneider 2004). Most of these detections were accomplished by 
the California-Carnegie group (Marcy, Butler, Fischer, & Vogt 2003). The iodine technique 
- thanks to its conceptual simplicity - - is the most commonly adopted way to obtain 
precision radial velocities. Iodine absorption cells are available on many spectrographs - 
HIRES at the 10m Keck I (Keck Observatory), Hamilton at the 3m Shane (Lick Obser- 
vatory), UCLES at the 3.9m Anglo- Australian Telescope (Anglo- Australian Observatory), 
HRS at the 9m HET (McDonald Observatory), MIKE at the 6.5m Magellan (Las Campanas 
Observatory), UVES at the 8m Kueyen (Cerro Paranal), HDS at the 8.2m Subaru (National 
Astronomical Observatory of Japan) and many others — and are used for planet detections. 
Unfortunately, in its classical version the iodine technique can only be applied to single 
stars or spectroscopic binaries whose secondaries are so faint that their spectral lines are 
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undetectable. This is dictated by the need to supply an observed template spectrum of each 
binary component for Eq. (1), but in the case of double-lined spectroscopic binaries it cannot 
be accomplished since their spectra are obviously always composite and time variable. 

In this paper, we present a novel technique which employs an I 2 absorption cell to obtain 
precision radial velocities of both components of double-lined spectroscopic binaries. The 
technique and its tests are described in Section 2 and its applications to stellar astronomy 
and searches for planets in binary and multiple stellar systems are discussed in Section 3. 

2. The Method 

We can measure precise radial velocities of both components of a spectroscopic binary 
with an I 2 absorption cell in the following way. First, we always take two subsequent expo- 
sures of each (binary) target — one with and the other without the I 2 cell. This is contrary 
to the standard approach for single stars where an exposure without the cell (a template) 
is taken only once. This way we obtain an instantaneous template which is used to model 
only the adjacent exposure taken with the cell. Next, we perform the usual least-squares fit 
and obtain the parameters described in Eq. (1). Obviously, the derived Doppler shift, AA^ 
(where % denotes the epoch of the observation), carries no meaning since each time a differ- 
ent template is used (besides it describes a Doppler "shift" of a composite spectrum that is 
typically different at each epoch). However, the parameters (in particular the wavelength 
solution and the parameters describing PSF) are accurately determined and can be used to 
extract the star spectrum, I*^ S (X), for each epoch i, by inverting the equation (1) 

' t(A) = [P obs W 0" 1 PSF'VT l2 (X) (2) 

where C§> _1 denotes deconvolution (carried out using a modified Jansson technique; Gilliland 
et al. 1992) and PSF 1 , symbolically, the set of parameters describing PSF at the epoch i. 
Such a star spectrum has an accurate wavelength solution, is free of the I 2 lines and the 
influence of a varying PSF. In the final step, the velocities of both components of a binary 
target can be measured with the well known two-dimensional cross-correlation technique 
TODCOR (Zucker & Mazeh 1994) using as templates the synthetic spectra derived with the 
ATLAS 9 and ATLAS 12 programs (Kurucz 1995) and matched to the observed spectrum, 
(Konacki, Torres, Sasselov, & Jha 2003a) I s (\). The formal errors of the velocities can be 
derived from the scatter between the velocities from different echelle orders or using the 
formalism of TODCOR (Zucker & Mazeh 1994). 

The above approach will produce radial velocities that may suffer from two possible 
sources of systematic errors. (1) The numerical procedure of extracting a star spectrum 
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from an exposure taken with the I2 cell could imaginably introduce systematic errors. (2) 
TODCOR itself could also produce systematic errors. It has never been tested at a high 
velocity precision (let us say, ~10 ms" 1 ) and the use of synthetic spectra could be a source 
of systematic errors. 

In order to demonstrate our method and investigate the systematic errors, we have 
observed a number of targets with HIRES spectrograph at the Keck I telescope. HIRES 
equipped with an I 2 absorption cell has been successfully used to detect many extrasolar 
planets (Marcy, Butler, Fischer, & Vogt 2003). We have observed HD 209458 (a single star 
with a known short period extrasolar planet, V=7.6 mag), HD 109358 (a velocity standard, 
V=4.3 mag), HD 4676 and HD 282975 (double-lined spectroscopic binaries in Pleiades, 
V=5.0 and 10.0 mag, respectively) and HD 206901 (a triple system, V=4.2 mag). The data 
span is approximately two years for HD 209458 (from mid 2002 to mid 2004), one year for 
HD 4676, HD 282975 and HD 206901 (from mid 2003 to mid 2004) and ~10 minutes for 
HD 109358 (on 29 Mar 2004). With the exception of HD 209458 whose spectra had a typical 
SNR of 40-80 (HD 209458 was initially used as a standard star for a separate transiting planet 
search and was observed with moderate exposure times; see Konacki, Torres, Sasselov, & Jha 
2003a; Konacki, Torres, Jha, & Sasselov 2003b), the SNR of all other spectra was typically 
in the range of 150-300. 

Since in our approach we do not assume that an observed object must have a composite 
spectrum, the method can be applied to single stars as well. Obviously, there is no need to 
analyze the spectra of single stars in such a way but the velocities computed in both ways 
(Eqs 1 and 2) can be compared and used to explore possible additional errors originating in 
the proposed procedure. Two such comparisons are presented in Figures 1 and 2. The first 
figure shows 6 velocity measurements of HD 109358 over the period of about 10 minutes. 
HD 109358 is a high-precision radial velocity standard proposed by the Geneva group (Udry, 
Mayor, & Queloz 1999). Clearly, the two velocity sets have very similar scatter and no 
additional errors at the 3 ms" 1 can be seen. The second test (Fig. 2) concerns HD 209458 
(the data span is two years), a well known star harboring a transiting giant planet with 
the orbital period of 3.52 days. Its orbital parameters are very accurately known (Naef et 
al. 2004) and hence were fixed in this test. The only parameter for which we fit was a 
velocity offset. The test does not reveal any additional errors either. In fact, the velocities 
computed using a cross-correlation function (CCF) have a smaller rms (10 vs 14 ms -1 ) than 
the velocities computed with the traditional iodine approach — most likely due to a mediocre 
SNR of the spectra of HD 209458 (the iodine technique with PSF modeling works best on 
high SNR spectra). 

In order to test our method on a truly binary spectrum we have observed HD 4676 



- 5 - 



(64 Psc) - - a nearby (23 pc) binary star composed of two F8 dwarfs. HD 4676AB has 
the orbital period of 13.8 days and the brightness ratio of ~ 0.9. Abt & Levy (1976) were 
the first to determine its spectroscopic orbit which was later improved by Nadal, Ginestet, 
Carquillat, & Pedoussaut (1979) and Duquennoy & Mayor (1991). Moreover, its astrometric 
orbit is known from the Palomar Testbed Interferometer (PTI; Boden et al. 1999). This 
has important consequences for our work as such astrometric data can be used to verify our 
precision radial velocities — orbital models for astrometric and spectroscopic measurements 
share most of the orbital parameters and can be used together to perform a combined least- 
squares fit (see e.g. Konacki & Lane 2004). 

We have observed HD 4676 over the period of one year on five occasions (11 Aug 2003, 12 
Aug 2003, 17 Nov 2003, 18 Nov 2003 and 18 Jul 2004) and obtained 24 radial velocity (RV) 
measurements for each component. The average formal RV error is 25 ms _1 for the primary 
and 28 ms _1 for the secondary. The quality of the data was verified by performing a combined 
fit to our RVs and the astrometric data published by Boden et al. (1999). The details of the 
fit are presented in Figure 3. The post-fit rms of 24 ms _1 is in a very reasonable agreement 
with the formal RV errors and is 27 times smaller than the rms reported by Duquennoy & 
Mayor (1991). It obviously results in a much better determination of the orbital and physical 
parameters of the system. We discuss it in detail in Section 3.1. 

The most important remaining issue is whether TODCOR and the use of synthetic 
spectra could introduce any systematic errors to RVs. While errors that simply increase 
the post-fit rms are not particularly troubling, one can in principle expect errors that may 
depend on the orbital phase and possibly affect the best-fit values of the orbital parameters. 
A common sense way to detect such errors is to apply our method to simulated spectra of 
HD 4676 and try to recover the RVs assumed in the simulation. We have performed such tests 
and have not detected any systematic errors. There is however an even better way to control 
these errors. As already mentioned, thanks to a high precision astrometry of HD 4676 from 
PTI (Boden et al. 1999), one can perform a combined fit to the astrometric and spectroscopic 
data. If there existed any significant phase-dependent systematic variations of the velocities, 
they would be detectable in the post-fit residuals from the combined solution. Clearly, these 
are absent in our solution (Fig. 3) and \ 2 1 DOF is 1.18. 

It should be noted that the best-fit rms of 24 ms _1 for the RVs of the primary and 
secondary is clearly larger than one would expect if the RV errors were photon-limited (2-3 
ms" 1 ; Butler et al. 1996). Possible sources of additional errors include: a mismatch between 
the synthetic templates and observed spectra and astrophysical phenomena (e.g. starspots, 
convective inhomogeneities, see Saar, Butler, & Marcy 1998). We do not have yet enough 
data for other binaries to fully address this issue. Nevertheless, it is reasonable to expect 



- 6- 



that for later type stars (with more spectral lines), the velocity precision will be higher. This 
appears to be confirmed by our velocities of HD 282975, a double-lined spectroscopic binary 
from Pleiades (P or b = 26 d ), which is composed of two G6 dwarfs. The post-fit rms is 11 
ms _1 for the primary and 13 ms" 1 for the secondary and is ~ 80 times smaller than the rms 
reported by Mermilliod, Rosvick, Duquennoy, & Mayor (1992). The best fit is based on only 
8 RVs but from our experience with HD 4676 we know that early accuracy is later confirmed 
by a larger data set. 

Finally, let us note that the iodine technique for single stars typically produces RVs 
which have an arbitrary offset (unless special steps are undertaken, see Nidever et al. 2002). 
In our approach which employs synthetic spectra one can in principle calibrate and refer the 
velocities to a standard system of stars. We plan to do it once we collect enough observations 
of the RV standards. 



3. Applications 

3.1. Stellar astronomy 

Modern tests of stellar structure and evolution models require stellar masses accurate to 
2% or better (Andersen 1991). Stellar masses are traditionally derived from the observations 
of binary stars by combining (1) photometric (light curves) and spectroscopic data (RVs) 
for double-lined eclipsing binaries, (2) relative astrometry and spectroscopic data for double- 
lined binaries, (3) from spectroscopy, relative astrometry and parallax for single-lined binaries 
and (4) absolute astrometry. There are currently known about 60 binaries with masses of 
the components accurate at the level of 2% or better (see e.g. Lastennet & Valls-Gabaud 
2002). Almost all of them are eclipsing binaries and most of them have components with 
the masses in the 1-3 M range (more massive and hence larger stars are more likely to be 
eclipsing, for a given orbital period). 

Optical interferometry offers the unique opportunity to determine visual orbits of double- 
lined spectroscopic binaries (apparent semi-major axes at the milliarcsecond level; for a recent 
discussion see Konacki & Lane 2004) and subsequently masses of the components, regardless 
of their size. Up to date, about twenty spectroscopic binaries have had their orbits deter- 
mined with optical interferometers (Quirrenbach 2001). In a few cases, when combined with 
radial velocity measurements, they already produced component masses accurate at the level 
of 1-2%. 

Our RVs constitute an excellent complement to precise astrometric orbits from optical 
interferometers. This is demonstrated well by HD 4676. Boden et al. (1999) combined 
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their interferometric data of HD 4676 from PTI with the radial velocity measurements of 
Duquennoy & Mayor (1991). They have determined the masses to be Mi = 1.223 ±O.O21M 
(1.7%), M 2 = 1.170 ± O.O18M (1.5%), the velocity amplitudes K x = 57.35 ± 0.31 kms" 1 , 
K 2 = 59.95 ±0.32 kms -1 and the inclination, i = 73.80 ±0.92 deg. The relative error in the 
masses (AM/M) is dominated by the error in the orbital inclination, 3Aicos(i)/ sin(i) 
1.4%, but the errors in the velocity amplitudes also contribute. In order to make the best 
out of PTI data, we used our RVs and performed a new combined fit. The best-fit orbital 
and physical parameters of the system are presented in Tables 1 and 2. The relative error in 
the masses is now 1.2% and the overall quality of the fit has improved. In this specific case, 
the ultimate accuracy of the masses is still limited by the accuracy of the orbital inclination. 
However, since the radial velocities are so precise, one can simply take more astrometric 
data with PTI and improve the masses even more. Clearly, our approach will be of great 
use for all those binaries that have orbital inclinations closer to an edge-on configuration. 
In particular, light curves of eclipsing binaries combined with precise RVs should produce 
masses accurate at the 0.1% level on a regular basis. 

Precise RVs can obviously be used not only to determine the masses of the companions 
but the distance to the binary as well. For example, the parallax of HD4676 is formally 
very accurate, 43.496 ± 0.089 mas (the distance of 22.991 ± 0.047 pc; the final accuracy will 
depend on how well the brightness ratio of the binary components is known). In this context, 
HD 282975 (Fig. 4) appears as a particularly interesting binary as it can be used to derive 
a very precise distance to Pleiades. To this end, one needs interferometric-astrometry of 
HD 282975 which should be easy to obtain with the Keck Interferometer (HD 282975 is too 
faint for PTI). Finally, precise RVs can be used to explore the regime of circular or almost 
circular orbits (the error of the eccentricity of HD 4676 is only 6.0 x 10~ 4 , Tab. 1) and the 
cut-off period below which orbits are circularized. 



3.2. Extrasolar planets in binary stellar systems 

Duquennoy & Mayor (1991) have demonstrated that the frequency of binaries (BF) 
among field stars older than 1 Gyr is 57%. The studies of multiplicity of pre main-sequence 
stars (PMS) in the Taurus and Ophiuchus star forming regions show that BF for systems in 
the separation range 1 to 150 AU is twice as large as that of the older field stars (Simon et al. 
1995). Further investigations have concluded that BF is lower for young stellar clusters (and 
similar to BF of the field stars) and that the binary frequency for PMS seems to be anti- 
correlated with the stellar density (Mathieu, Ghez, Jensen, & Simon 2000). Nonetheless, 
BF is very high for both field and pre main-sequence stars and one can argue that it may 
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not be possible to assess the overall frequency and properties of extrasolar planets without 
addressing binary (and multiple) stellar systems. 

Supporting arguments come from the presence of circumstellar and circumbinary disks 
around binaries. One of the prime examples of the circumstellar disks in a binary system is 
the case of LI 551 IRS 5. Rodriguez et al. (1998) show that LI 551 IRS 5 is a binary PMS 
with the separation of 45 AU in which each component is surrounded by a disk. The radii of 
the disks are 10 AU and the estimated masses are 0.06 and 0.03 M Q , supposedly enough to 
produce planets. Recently, McCabe, Duchene, & Ghez (2003) have spatially resolved mid- 
infrared scattered light from the protoplanetary disk around the secondary of the PMS binary 
HK Tau AB. The inferred sizes of the dust grains are in the range 1.5-3 /zm which suggests 
that the first step in the planet formation, the dust grain growth, has occurred in this disk. 
Millimeter and submillimeter measurements of the dust continuum emission enable us to 
measure the total masses of disks. These observations show that the circumbinary disks may 
be reduced in size and mass but still are present even in close systems. The circumbinary 
disks are observed at millimeter wavelengths around many PMS spectroscopic binaries. Such 
massive disks are however rare around wide binaries with separations of 1-100 AU. This is 
confirmed by numerical works that predict circumstellar and circumbinary disks truncated 
by the companions (Lubow & Artymowicz 2000). The circumstellar disks have the outer 
radii 0.2-0.5 times the binary separation while the circumbinary disks have the inner radii 
2-3 times the semi-major axis of the binary. Finally, the measurements of the infrared excess 
emission show no difference in frequency of the excess among binaries and single stars. It 
indicates that the circumstellar material in binary systems may be similar in temperature 
and surface density to that in the disks surrounding single stars (Mathieu, Ghez, Jensen, & 
Simon 2000). 

The problem of stability of the planetary orbits in binaries was recognized a long time 
ago. It was mostly approached with the numerical studies of the elliptic restricted three-body 
problem. The orbital configurations considered include the so-called P-type (Planet-type, 
circumbinary orbits), S-type (Satellite-type, circumprimary or circumsecondary orbits) and 
L-type orbits (Librator-type, orbits around stable Lagrangian points L4 or L5 for the mass 
ratios fi < 0.04). There are many papers concerning the stability of S-type motions (e.g. 
Benest 2003; Pilat-Lohinger & Dvorak 2002; Rabl & Dvorak 1988). They aimed at developing 
empirical stability criteria in the framework of the circular three-body problem (see e.g. 
Graziani & Black 1981; Black 1982; Pendleton & Black 1983), the stability of periodic orbits 
and the stability of the test particles in binary systems as a function of the eccentricity of 
the binary. The P-type motions were also investigated (Pilat-Lohinger, Funk, & Dvorak 
2003; Broucke 2001; Holman & Wiegert 1999); as were the L-type orbits (see e.g. Laughlin 
& Chambers 2002). Unfortunately, most of the analytical papers deal only with circular 
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binary orbits, numerical studies are confined to special mass ratios and the integration times 
used are relatively short. Also, they are almost exclusively restricted to the framework of 
the three body problem. Some of these issues are addressed by Holman & Wiegert (1999) 
who studied a range of mass ratios, eccentricities and used long integration times (at least 
10 4 periods of the binary). Nonetheless, we can safely conclude that planets in binary stellar 
systems, once formed, would enjoy a wide range of stable orbits. 

However, the theories of planet formation in binary stellar systems are still at early 
stages. Whitmire, Matese, Criswell, & Mikkola (1998) studied terrestrial planet growth in 
the circumprimary habitable zones in a binary system. They considered a 4-body system of 2 
stars and 2 planetesimals for which by varying the orbital parameters of the binary (the semi- 
major axis, eccentricity and mass ratio) they were able to determine a critical semi-major 
axis of the binary below which the secondary does not allow a growth of the planetesimals 
(planetesimals are accelerated by the secondary, their relative velocity is larger than criti- 
cal and their collisions become destructive). Based on this criterion, they concluded that 
about 60% of nearby solar-type binaries cannot be excluded from having a habitable planet. 
Marzari & Scholl (2000) analyzed a Cen AB (the semi-major axis of 23 AU, eccentricity 
of 0.52, mass ratio 1.1/0.92), a prototype close binary system, and demonstrated that the 
planetesimals can accrete into planetary embryos. Barbieri, Marzari, & Scholl (2002) contin- 
ued the study and showed that planetary embryos can grow into terrestrial planets in about 
50 Myr. Somewhat contrary to this result, Nelson (2000), who analyzed a binary system 
similar to L1551 IRS 5, found that the planet formation is unlikely in equal mass binary 
systems with the separation of about 50 AU. Yet another result by Boss (1998) claims that a 
stellar companion can induce a multi- Jupiter-mass planet formation. Clearly, there is a lack 
of consensus and the planet formation theories would certainly benefit from observational 
constraints. 

Unfortunately, current RV searches for extrasolar planets have been limited to single 
stars or binaries with large apparent separations. Despite that prejudice out of ~130 ex- 
trasolar planets, eighteen belong to stars which have also stellar companions (Eggenberger, 
Udry, & Mayor 2004). They even seem to posses characteristics different from those orbit- 
ing single stars. As recently pointed out by Eggenberger, Udry, & Mayor (2004), (1) the 
most massive short period planets orbit stars from multiple stellar systems and (2) all these 
planets (P or b < 40 days) tend to have very low eccentricities. It is a surprisingly interesting 
result given that current RV surveys suffer from a selection effect — the binary stars with 
separations smaller than ~2 arcsec have been excluded from the surveys. The exception 
here is HD 41004AB which has a separation of 0.5 arcsec. Presumably, this star entered 
the Geneva RV program as a single star and was later found to be a triple (sic!) system in 
which the primary has apparently a giant planet (Eggenberger, Udry, & Mayor 2004) and 
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the secondary and tertiary are very faint (an M star and a putative brown dwarf). Thus 
these stars from binary or multiple stellar systems which we know that harbor planets have 
a distant or/and a faint stellar companion. There is a need to perform an unbiased survey 
for extrasolar planets and substellar companions to stars from multiple systems which will 
probe different physical separations and brightness ratios. 

To this end, two years ago we have initiated a survey to detect planets in binary stellar 
systems with the Keck I and HET telescopes. Our survey covers about ~450 binary stars 
from the northern and southern hemispheres. Preliminary results for one of our represen- 
tative targets, re Peg are shown in Fig. 5. re Peg (HR 8315, HD 206901, HIP 107354) is a 
known triple system of which the spectral lines of the two brightest components AB (appar- 
ent semi-major axis 0.4 arcsec) can be easily identified (Am = 0.04, P orb = 11.59 yr) and 
the component B is itself a single-lined spectroscopic binary with the orbital period of ~5.97 
days (e.g. Mayor & Mazeh 1987). The radial velocities of the primary do not reveal any 
interesting variations and are of inferior quality due to a large V Sill % ~ 30 kms 1 . On the 
other hand, the RVs of the sharp-lined secondary are more than acceptable (given an early 
spectral type of the star, F3IV) with the rms of ~ 30 ms -1 . It is an improvement of 30 times 
in the rms compared to the work of Mayor & Mazeh (1987) and demonstrates the usefulness 
of our approach. Since the iodine cells are available on many high-resolution spectrographs, 
our method clearly opens exciting opportunities in the studies of other multiple systems. 

4. Summary 

We have introduced a new method to measure precision radial velocities of double-lined 
spectroscopic binaries with an iodine absorption cell. Our initial data sets demonstrate 20-30 
ms -1 velocity precision for "early" (F3-F8) type and ~10 ms -1 for later type binaries. If 
combined with interferometric astrometry, such velocities enable us to determine the masses 
of the components of binary stars accurate at least at the level of 1%. Our method also 
makes it possible to perform a search for extrasolar planets in binary and multiple stellar 
systems where the brightness of the primary and the secondary is comparable. This allows 
us to probe a regime of the parameter space of the formation and evolution of extrasolar 
planets not covered by other surveys. 
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Fig. 1. — Radial velocities of HD 109358 computed with (a) the traditional iodine (I 2 ) 
technique and with (6) the proposed method. 

Fig. 2. — Radial velocities of HD 209458 computed with (a) the traditional iodine (I 2 ) 
technique and with (b) the proposed method. The data span is two years. (c,d) Best-fit 
residuals as a function of phase. 

Fig. 3. — Observed (filled circles for the primary and open circles for the secondary) and 
modeled (solid line) radial velocities of HD 4676. (b,c) Best-fit residuals as a function of 
phase and time, (d) Histogram of the residuals. 

Fig. 4. — Observed (filled circles for the primary and open circles for the secondary) and 
modeled (solid line) radial velocities of HD 282975. (b,c) Best-fit residuals as a function of 
phase and time. 

Fig. 5. — Observed (filled circles) and modeled (solid line) radial velocities of HD 206901B. 
(b,c) Best-fit residuals as a function of phase and time, (d) Histogram of the residuals. 
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Table 1. Best-fit Orbital Parameters for HD 4676^ 



Parameter HD 4676 



Apparent semi-major axis, a (mas) 6.545 ± 0.0133 

Period, P (d) 13.8244906 ± 4.3xl0~ 5 

Time of periastron, T p (MJD) 50905.9746 ± 0.0067 

Eccentricity, e 0.23657 ± 0.00063 

Longitude of the periastron, u (deg) 203.057 ± 0.073 

Longitude of the ascending node, Q (deg) 207.41 ± 0.65 

Inclination, % (deg) 73.92 ± 0.80 

Magnitude difference (K band, assumed), AK . 0.11 

Velocity amplitude of the primary, K\ (km/s) . 57.552 ± 0.037 

Velocity amplitude of the secondary, K 2 (km/s) 59.557 ± 0.038 

Reduced x 2 , X 2 /DOF 1.18 



a Figures in parentheses are the la formal errors 



Table 2. Physical Parameters for HD 4676 



Parameter Primary Secondary 



Semi-major axis, ai, 2 (AU) 0.073953 ± 4.8 xlO" 5 0.076529 ± 5.0xl0~ 5 

Mass, M (M ) 1.210 ± 0.014 1.169 ± 0.014 

Parallax, k (mas) 43.496 ± 0.089 

Distance, d (pc) 22.991 ± 0.047 

Spectral type F8V F8V 



a Figures in parentheses are the la formal errors 



